BSXplorer: analytical framework for exploratory analysis of BS-seq data.

BACKGROUND: Bisulfite sequencing detects and quantifies DNA methylation patterns, contributing to our understanding of gene expression regulation, genome stability maintenance, conservation of epigenetic mechanisms across divergent taxa, epigenetic inheritance and, eventually, phenotypic variation. Graphical representation of methylation data is crucial in exploring epigenetic regulation on a genome-wide scale in both plants and animals. This is especially relevant for non-model organisms with poorly annotated genomes and/or organisms where genome sequences are not yet assembled on chromosome level. Despite being a technology of choice to profile DNA methylation for many years now there are surprisingly few lightweight and robust standalone tools available for efficient graphical analysis of data in non-model systems. This significantly limits evolutionary studies and agrigenomics research. BSXplorer is a tool specifically developed to fill this gap and assist researchers in explorative data analysis and in visualising and interpreting bisulfite sequencing data more easily. RESULTS: BSXplorer provides in-depth graphical analysis of sequencing data encompassing (a) profiling of methylation levels in metagenes or in user-defined regions using line plots and heatmaps, generation of summary statistics charts, (b) enabling comparative analyses of methylation patterns across experimental samples, methylation contexts and species, and (c) identification of modules sharing similar methylation signatures at functional genomic elements. The tool processes methylation data quickly and offers API and CLI capabilities, along with the ability to create high-quality figures suitable for publication. CONCLUSIONS: BSXplorer facilitates efficient methylation data mining, contrasting and visualization, making it an easy-to-use package that is highly useful for epigenetic research.

A Combined Effect of Polygenic Scores and Environmental Factors on Individual Differences in Depression Level.

The risk of depression could be evaluated through its multifactorial nature using the polygenic score (PGS) approach. Assuming a "clinical continuum" hypothesis of mental diseases, a preliminary assessment of individuals with elevated risk for developing depression in a non-clinical group is of high relevance. In turn, epidemiological studies suggest including social/lifestyle factors together with PGS to address the "missing heritability" problem. We designed regression models, which included PGS using 27 SNPs and social/lifestyle factors to explain individual differences in depression levels in high-education students from the Volga-Ural region (VUR) of Eurasia. Since issues related to population stratification in PGS scores may lead to imprecise variant effect estimates, we aimed to examine a sensitivity of PGS calculated on summary statistics of depression and neuroticism GWAS from Western Europeans to assess individual proneness to depression levels in the examined sample of Eastern Europeans. A depression score was assessed using the revised version of the Beck Depression Inventory (BDI) in 1065 young adults (age 18-25 years, 79% women, Eastern European ancestry). The models based on weighted PGS demonstrated higher sensitivity to evaluate depression level in the full dataset, explaining up to 2.4% of the variance (p = 3.42 × 10-7); the addition of social parameters enhanced the strength of the model (adjusted r2 = 15%, p < 2.2 × 10-16). A higher effect was observed in models based on weighted PGS in the women group, explaining up to 3.9% (p = 6.03 × 10-9) of variance in depression level assuming a combined SNPs effect and 17% (p < 2.2 × 10-16)-with the addition of social factors in the model. We failed to estimate BDI-measured depression based on summary statistics from Western Europeans GWAS of clinical depression. Although regression models based on PGS from neuroticism (depression-related trait) GWAS in Europeans were associated with a depression level in our sample (adjusted r2 = 0.43%, p = 0.019-for unweighted model), the effect was mainly attributed to the inclusion of social/lifestyle factors as predictors in these models (adjusted r2 = 15%, p < 2.2 × 10-16-for unweighted model). In conclusion, constructed PGS models contribute to a proportion of interindividual variability in BDI-measured depression in high-education students, especially women, from the VUR of Eurasia. External factors, including the specificity of rearing in childhood, used as predictors, improve the predictive ability of these models. Implementation of ethnicity-specific effect estimates in such modeling is important for individual risk assessment.

Pygenomics: manipulating genomic intervals and data files in Python.

SUMMARY: We present pygenomics, a Python package for working with genomic intervals and bioinformatic data files. The package implements interval operations, provides both API and CLI, and supports reading and writing data in widely used bioinformatic formats, including BAM, BED, GFF3, and VCF. The source code of pygenomics is provided with in-source documentation and type annotations and adheres to the functional programming paradigm. These features facilitate seamless integration of pygenomics routines into scripts and pipelines. The package is implemented in pure Python using its standard library only and contains the property-based testing framework. Comparison of pygenomics with other Python bioinformatic packages with relation to features and performance is presented. The performance comparison covers operations with genomic intervals, read alignments, and genomic variants and demonstrates that pygenomics is suitable for computationally effective analysis. AVAILABILITY AND IMPLEMENTATION: The source code is available at https://gitlab.com/gtamazian/pygenomics.

A collection of read depth profiles at structural variant breakpoints.

SWaveform, a newly created open genome-wide resource for read depth signal in the vicinity of structural variant (SV) breakpoints, aims to boost development of computational tools and algorithms for discovery of genomic rearrangement events from sequencing data. SVs are a dominant force shaping genomes and substantially contributing to genetic diversity. Still, there are challenges in reliable and efficient genotyping of SVs from whole genome sequencing data, thus delaying translation into clinical applications and wasting valuable resources. SWaveform includes a database containing ~7 M of read depth profiles at SV breakpoints extracted from 911 sequencing samples generated by the Human Genome Diversity Project, generalised patterns of the signal at breakpoints, an interface for navigation and download, as well as a toolbox for local deployment with user's data. The dataset can be of immense value to bioinformatics and engineering communities as it empowers smooth application of intelligent signal processing and machine learning techniques for discovery of genomic rearrangement events and thus opens the floodgates for development of innovative algorithms and software.

Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer.

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Long non-coding RNA-derived peptides are immunogenic and drive a potent anti-tumour response.

Protein arginine methyltransferase (PRMT) 5 is over-expressed in a variety of cancers and the master transcription regulator E2F1 is an important methylation target. We have explored the role of PRMT5 and E2F1 in regulating the non-coding genome and report here a striking effect on long non-coding (lnc) RNA gene expression. Moreover, many MHC class I protein-associated peptides were derived from small open reading frames in the lncRNA genes. Pharmacological inhibition of PRMT5 or adjusting E2F1 levels qualitatively altered the repertoire of lncRNA-derived peptide antigens displayed by tumour cells. When presented to the immune system as either ex vivo-loaded dendritic cells or expressed from a viral vector, lncRNA-derived peptides drove a potent antigen-specific CD8 T lymphocyte response, which translated into a significant delay in tumour growth. Thus, lncRNA genes encode immunogenic peptides that can be deployed as a cancer vaccine.

Host Genetic Variants Linked to COVID-19 Neurological Complications and Susceptibility in Young Adults-A Preliminary Analysis.

To date, multiple efforts have been made to use genome-wide association studies (GWAS) to untangle the genetic basis for SARS-CoV-2 infection susceptibility and severe COVID-19. However, data on the genetic-related effects of SARS-CoV-2 infection on the presence of accompanying and long-term post-COVID-19 neurological symptoms in younger individuals remain absent. We aimed to examine the possible association between SNPs found in a GWAS of COVID-19 outcomes and three phenotypes: SARS-CoV-2 infection, neurological complications during disease progression, and long-term neurological complications in young adults with a mild-to-moderate disease course. University students (N = 336, age 18-25 years, European ancestry) with or without COVID-19 and neurological symptoms in anamnesis comprised the study sample. Logistic regression was performed with COVID-19-related phenotypes as outcomes, and the top 25 SNPs from GWAS meta-analyses and an MR study linking COVID-19 and cognitive deficits were found. We replicated previously reported associations of the FURIN and SLC6A20 gene variants (OR = 2.36, 95% CI 1.31-4.24) and OR = 1.94, 95% CI 1.08-3.49, respectively) and remaining neurological complications (OR = 2.12, 95% CI 1.10-4.35 for SLC6A20), while NR1H2 (OR = 2.99, 95% CI 1.39-6.69) and TMPRSS2 (OR = 2.03, 95% CI 1.19-3.50) SNPs were associated with neurological symptoms accompanying COVID-19. Our findings indicate that genetic variants related to a severe COVID-19 course in adults may contribute to the occurrence of neurological repercussions in individuals at a young age.

A Genome-Wide Association Study Reveals a BDNF-Centered Molecular Network Associated with Alcohol Dependence and Related Clinical Measures.

At least 50% of factors predisposing to alcohol dependence (AD) are genetic and women affected with this disorder present with more psychiatric comorbidities, probably indicating different genetic factors involved. We aimed to run a genome-wide association study (GWAS) followed by a bioinformatic functional annotation of associated genomic regions in patients with AD and eight related clinical measures. A genome-wide significant association of rs220677 with AD (p-value = 1.33 × 10-8 calculated with the Yates-corrected χ2 test under the assumption of dominant inheritance) was discovered in female patients. Associations of AD and related clinical measures with seven other single nucleotide polymorphisms listed in previous GWASs of psychiatric and addiction traits were differently replicated in male and female patients. The bioinformatic analysis showed that regulatory elements in the eight associated linkage disequilibrium blocks define the expression of 80 protein-coding genes. Nearly 68% of these and of 120 previously published coding genes associated with alcohol phenotypes directly interact in a single network, where BDNF is the most significant hub gene. This study indicates that several genes behind the pathogenesis of AD are different in male and female patients, but implicated molecular mechanisms are functionally connected. The study also reveals a central role of BDNF in the pathogenesis of AD.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration.

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

An account of Fusarium wilt resistance in flax Linum usitatissimum: The disease severity data.

A collection of flax accessions maintained by the Russian Federal Research Center for Bast Fiber Crops was characterized to evaluate its resistance to Fusarium wilt. 297 samples representing different morphotypes and selection status were infected with highly virulent MI39 strain of Fusarium oxysporum f. sp. lini. Evaluation of disease symptoms was performed at the full emergence stage in a greenhouse. The experiment lasted for 3 successive years. The disease severity index (DSI) record was obtained for every genotype in each year of the experiment. The data set was produced in a framework of a project focused on both a) deciphering the mechanisms of plant immunity, and b) development of flax cultivars resistant to Fusarium wilt. The data is available via Figshare repository.

The Genetic Landscape of Fiber Flax.

Genetic diversity in a breeding program is essential to overcome modern-day environmental challenges faced by humanity and produce robust, resilient crop cultivars with improved agronomic characteristics, as well as to trace crop domestication history. Flax (Linum usitatissimum), one of the first crops domesticated by mankind, has been traditionally cultivated for fiber as well as for medicinal purposes and as a nutritional product. The origins of fiber flax are hidden in the mists of time and can be hypothetically traced back to either the Indo-Afghan region or Fertile Crescent. To shed new light on fiber flax genetic diversity and breeding history, in this study, we presented a comprehensive analysis of the core collection of flax (306 accessions) of different morphotypes and geographic origins maintained by the Russian Federal Research Center for Bast Fiber Crops. We observed significant population differentiation between oilseed and fiber morphotypes, as well as mapped genomic regions affected by recent breeding efforts. We also sought to unravel the origins of kryazhs, Russian heritage landraces, and their genetic relatedness to modern fiber flax cultivars. For the first time, our results provide strong genetic evidence in favor of the hypothesis on kryazh's mixed origin from both the Indo-Afghan diversity center and Fertile Crescent. Finally, we showed predominant contribution from Russian landraces and kryazhs into the ancestry of modern fiber flax varieties. Taken together, these findings may have practical implications on the development of new improved flax varieties with desirable traits that give farmers greater choice in crop management and meet the aspirations of breeders.

Stochastic Effects in Retrotransposon Dynamics Revealed by Modeling under Competition for Cellular Resources.

Transposons are genomic elements that can relocate within a host genome using a 'cut'- or 'copy-and-paste' mechanism. They make up a significant part of many genomes, serve as a driving force for genome evolution, and are linked with Mendelian diseases and cancers. Interactions between two specific retrotransposon types, autonomous (e.g., LINE1/L1) and nonautonomous (e.g., Alu), may lead to fluctuations in the number of these transposons in the genome over multiple cell generations. We developed and examined a simple model of retrotransposon dynamics under conditions where transposon replication machinery competed for cellular resources: namely, free ribosomes and available energy (i.e., ATP molecules). Such competition is likely to occur in stress conditions that a malfunctioning cell may experience as a result of a malignant transformation. The modeling revealed that the number of actively replicating LINE1 and Alu elements in a cell decreases with the increasing competition for resources; however, stochastic effects interfere with this simple trend. We stochastically simulated the transposon dynamics in a cell population and showed that the population splits into pools with drastically different transposon behaviors. The early extinction of active Alu elements resulted in a larger number of LINE1 copies occurring in the first pool, as there was no competition between the two types of transposons in this pool. In the other pool, the competition process remained and the number of L1 copies was kept small. As the level of available resources reached a critical value, both types of dynamics demonstrated an increase in noise levels, and both the period and the amplitude of predator-prey oscillations rose in one of the cell pools. We hypothesized that the presented dynamical effects associated with the impact of the competition for cellular resources inflicted on the dynamics of retrotransposable elements could be used as a characteristic feature to assess a cell state, or to control the transposon activity.

Genomic Regions Associated with Fusarium Wilt Resistance in Flax.

Modern flax cultivars are susceptible to many diseases; arguably, the most economically damaging of these is the Fusarium wilt fungal disease. Over the past decades international flax breeding initiatives resulted in the development of resistant cultivars. However, much remains to be learned about the mechanisms of resistance to Fusarium infection in flax. As a first step to uncover the genetic factors associated with resistance to Fusarium wilt disease, we performed a genome-wide association study (GWAS) using 297 accessions from the collection of the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. These genotypes were infected with a highly pathogenic Fusarium oxysporum f.sp. lini MI39 strain; the wilt symptoms were documented in the course of three successive years. Six different single-locus models implemented in GAPIT3 R package were applied to a selected subset of 72,526 SNPs. A total of 15 QTNs (Quantitative Trait Nucleotides) were detected during at least two years of observation, while eight QTNs were found during all three years of the experiment. Of these, ten QTNs occupied a region of 640 Kb at the start of chromosome 1, while the remaining QTNs mapped to chromosomes 8, 11 and 13. All stable QTNs demonstrate a statistically significant allelic effect across 3 years of the experiment. Importantly, several QTNs spanned regions that harbored genes involved in the pathogen recognition and plant immunity response, including the KIP1-like protein (Lus10025717) and NBS-LRR protein (Lus10025852). Our results provide novel insights into the genetic architecture of flax resistance to Fusarium wilt and pinpoint potential candidate genes for further in-depth studies.

A comprehensive dataset of flax (Linum uitatissimum L.) phenotypes.

A collection of flax accessions from Russian Federal Research Center for Bast Fiber Crops was characterised to evaluate its phenotypic diversity. 406 samples representing different morphotypes were selected for thorough quantitative assessment of various agronomic traits. We measured height, length of technical part of the stem, technical part weight, inflorescence length, number of bolls and seeds per plant, 1000 seed weight, the diameter of the stem, the number of internodes and finally, distance between internodes. The fiber quality was estimated by calculating stem slenderness, stem taperingness and elementary fiber length. The dataset was produced in a framework of a project focused on characterization of diversity of flax genotypes and phenotypes, as well as on identification of genomic regions associated with various traits, it is hosted on Figshare.

A Genomic Blueprint of Flax Fungal Parasite Fusarium oxysporum f. sp. lini.

Fusarium wilt of flax is an aggressive disease caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. lini. It is a challenging pathogen presenting a constant threat to flax production industry worldwide. Previously, we reported chromosome-level assemblies of 5 highly pathogenic F. oxysporum f. sp. lini strains. We sought to characterize the genomic architecture of the fungus and outline evolutionary mechanisms shaping the pathogen genome. Here, we reveal the complex multi-compartmentalized genome organization and uncover its diverse evolutionary dynamics, which boosts genetic diversity and facilitates host adaptation. In addition, our results suggest that host of functions implicated in the life cycle of mobile genetic elements are main contributors to dissimilarity between proteomes of different Fusaria. Finally, our experiments demonstrate that mobile genetics elements are expressed in planta upon infection, alluding to their role in pathogenicity. On the whole, these results pave the way for further in-depth studies of evolutionary forces shaping the host-pathogen interaction.

A genome-wide association study identifies a gene network associated with paranoid schizophrenia and antipsychotics-induced tardive dyskinesia.

In the present study we conducted a genome-wide association study (GWAS) in a cohort of 505 patients with paranoid schizophrenia (SCZ), of which 95 had tardive dyskinesia (TD), and 503 healthy controls. Using data generated by the PsychENCODE Consortium (PEC) and other bioinformatic databases, we revealed a gene network, implicated in neurodevelopment and brain function, associated with both these disorders. Almost all these genes are in gene or isoform co-expression PEC network modules important for the functioning of the brain; the activity of these networks is also altered in SCZ, bipolar disorder and autism spectrum disorders. The associated PEC network modules are enriched for gene ontology terms relevant to the brain development and function (CNS development, neuron development, axon ensheathment, synapse, synaptic vesicle cycle, and signaling receptor activity) and to the immune system (inflammatory response). Results of the present study suggest that orofacial and limbtruncal types of TD seem to share the molecular network with SCZ. Paranoid SCZ and abnormal involuntary movements that indicate the orofacial type of TD are associated with the same genomic loci on chromosomes 3p22.2, 8q21.13, and 13q14.2. The limbtruncal type of TD is associated with a locus on chromosome 3p13 where the best functional candidate is FOXP1, a high-confidence SCZ gene. The results of this study shed light on common pathogenic mechanisms for SCZ and TD, and indicate that the pathogenesis of the orofacial and limbtruncal types of TD might be driven by interacting genes implicated in neurodevelopment.

Correction to: Characterising a human endogenous retrovirus(HERV)-derived tumour-associated antigen: enriched RNA-Seq analysis of HERV-K(HML-2) in mantle cell lymphoma cell lines.

An amendment to this paper has been published and can be accessed via the original article.

PRMT5 promotes cancer cell migration and invasion through the E2F pathway.

The pRb-E2F pathway is a critical point of regulation in the cell cycle and loss of control of the pathway is a hallmark of cancer. E2F1 is the major target through which pRb exerts its effects and arginine methylation by PRMT5 plays a key role in dictating E2F1 activity. Here we have explored the functional role of the PRMT5-E2F1 axis and highlight its influence on different aspects of cancer cell biology including viability, migration, invasion and adherence. Through a genome-wide expression analysis, we identified a distinct set of genes under the control of PRMT5 and E2F1, including some highly regulated genes, which influence cell migration, invasio and adherence through a PRMT5-dependent mechanism. Most significantly, a coincidence was apparent between the expression of PRMT5 and E2F1 in human tumours, and elevated levels of PRMT5 and E2F1 correlated with poor prognosis disease. Our results suggest a causal relationship between PRMT5 and E2F1 in driving the malignant phenotype and thereby highlight an important pathway for therapeutic intervention.

The Genome Sequence of Five Highly Pathogenic Isolates of Fusarium oxysporum f. sp. lini.

Fusarium wilt is the most destructive fungal disease in flax, limiting flax cultivation in all the main flax and linseed growing countries. The causative agent is seedborne and soilborne fungus F. oxysporum f. sp. lini. Here, we report, for the first time, genome assemblies of five highly pathogenic isolates of Fusarium oxysporum f. sp. lini, namely monoisolate 39 and strains F329, F324, F282, F287. In addition, syntenic analysis provided a powerful approach to distinguish between core and lineage-specific parts of the genome. These results lay a solid foundation for comparative genomics studies of plant fungal pathogens, evolution of pathogenicity, and virulence factors underlying the dynamics of host-pathogen interactions, thus eventually offering solutions to Fusarium disease control.

Current challenges and possible future developments in personalized psychiatry with an emphasis on psychotic disorders.

A personalized medicine approach seems to be particularly applicable to psychiatry. Indeed, considering mental illness as deregulation, unique to each patient, of molecular pathways, governing the development and functioning of the brain, seems to be the most justified way to understand and treat disorders of this medical category. In order to extract correct information about the implicated molecular pathways, data can be drawn from sampling phenotypic and genetic biomarkers and then analyzed by a machine learning algorithm. This review describes current difficulties in the field of personalized psychiatry and gives several examples of possibly actionable biomarkers of psychotic and other psychiatric disorders, including several examples of genetic studies relevant to personalized psychiatry. Most of these biomarkers are not yet ready to be introduced in clinical practice. In a next step, a perspective on the path personalized psychiatry may take in the future is given, paying particular attention to machine learning algorithms that can be used with the goal of handling multidimensional datasets.